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Using a model of force balance in Saturn's disk-like magnetosphere, we show 
that variations in hot plasma pressure can change the magnetic field con- 
figuration. This effect changes (i) the location of the magnetopause, even at 
fixed solar wind dynamic pressure, and (ii) the magnetic mapping between 
ionosphere and disk. The model uses equatorial observations as a boundary 
condition — we test its predictions over a wide latitude range by comparison 
with a Cassini high-inclination orbit of magnetic field and hot plasma pres- 
sure data. We find reasonable agreement over time scales larger than the pe- 
riod of Saturn kilometric radiation (also known as the camshaft period). 
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1. Introduction 

Saturn's equatorial, rotating plasma disk is threaded by a magnetic field which changes 
from a quasi-dipolar geometry, at distances <12R S (Saturn's radius Rs = 60280 km), 
to a "magnetodisk" in the outer region, where the radial field component is dominant 
at distances ~1-2R S from the equator (e.g Arridge et al. [2007]). This magnetodisk is 
supported by an azimuthal ring current whose solenoid-like field adds to the planetary 
dipole. In early studies, the ring current was represented by a bound annular region with 
four free parameters: a uniform half-thickness, inner and outer radii, and a current scaling 
parameter (e.g. Connerney et al. [1981]). In this model, current density was assumed 
inversely proportional to cylindrical radial distance p. This formalism ("CAN81 disk") 
was recently adapted to analyze Cassini magnetic data, establishing how ring current 
parameters depend on magnetopause size [Bunce et al, 2007]. 

More recently, a formalism for calculating self-consistent field and plasma distributions 
for the Jovian magnetodisk [Caudal, 1986] has been implemented for Saturn (Achilleos 
et al. [2010], hereafter A 10). This approach integrates equations of force balance through- 
out a rotating, axisymmetric magnetosphere (§2). In §3, we use the A 10 model to show 
that the variability of hot plasma pressure, as observed, affects the magnetospheric field 
structure — increased hot plasma content generally producing a more "inflated" or "disk- 
like" field geometry. To test the model's high-latitude structure, we compare its outputs 
with field and particle data from the Rev 40 high-inclination orbit of Cassini. We discuss 
the agreement between model and data and summarize our conclusions in §4. 

2. Magnetodisk Model Inputs 
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We have improved AlO's representation of the cold equatorial ion temperature over 
10-25 Rs, by combining temperatures for protons and water group ions from McAndrews 
et al. [2009] (-10-30 R s ) and Wilson et al. [2008] (5-10 R s ), who analyzed energy spectra 
acquired by the ion mass spectrometer (IMS), a subsystem of the Cassini Plasma Spec- 
trometer (CAPS) [Young et al, 2004]. Temperature moments versus equatorial distance, 
along with quartic fits, are shown in Fig. la,b (fit coefficients available on request). AWs 
initial study used only Wilson et al. [2008] 's temperatures, with parabolic extrapolation 
beyond 10 Rs- Using the updated profiles, we derived the single plasma temperature 
T c (p) required by the model, whose field lines are loci of constant, isotropic temperature. 
This was done as in A10 , by taking weighted average moments for each ion species, then 
combining these according to a prescribed plasma composition. 

To model hot plasma content, we used the equatorial product Ph(p) V a (p) of hot plasma 
pressure (assumed constant along a field line) and unit flux tube volume (Caudal [1986], 
A 10). We extended the Cassini Magnetospheric Imaging Instrument (MIMI) pressure 
dataset used in AlO's initial study [Sergis et al, 2007] to cover 23 near-equatorial orbits 
in total, between October, 2004 and March, 2006. The hot pressure arises from ions 
(mainly H+ and 0+) with energies >3keV (e.g. Sergis et al. [2009]). The moments 
are computed in the spacecraft frame, from an incomplete sampling of the phase space. 
However, for p>12Rg, energetic ions have typical thermal speeds —10 times higher than 
the plasma flow speed. The resulting hot pressures are uncertain by factors of —2—3, 
within the scatter of the presented data. We used pressure moments sampled every 10 
minutes, taken from the dayside, equatorial magnetosphere: (i) absolute latitudes <5°, 
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(ii) Saturn local times 09:00-15:00, and (iii) locations between each orbit's last inbound 
and next outbound magnetopause crossings. 

All hot pressures are shown as gray symbols in Fig. lc, collapsed onto a single profile in 
p. For comparison, data for Cassini Rev 3 (February 14-23, 2005) are overlain in black 
(data from this orbit were previously analyzed by A 10). Fig. Id shows the median, 10% 
and 90% quantiles. Solid curves are cubic fits to these parameters, which correspond to 
"average", "disturbed" and "quiescent" ring current states (following Sergis et al. [2007] 
and A 10). We assumed a zero-slope profile past ~23 R s (dashed lines in Fig. Id), however 
this did not affect our conclusions. Repeating this parametrization, incorporating all local 
times, elevated the quantile pressures by ~40%. 

3. Comparison of Magnetodisk Models and Observations 

Fig. 2 illustrates the model magnetodisk's response to ring current state. Left and 
right columns show calculations for fixed magnetopause radius Rmp=25 Rs (average value 
for Saturn [Achilleos et al, 2008]), but for quiescent and disturbed states (Fig. Id). 
Fig. 2a shows that field lines in the disturbed model become more radially "stretched" 
for p~10-15Rs. The labeled field lines intersect northern ionospheric colatitudes ^=15°, 
17° and 19°. These values are first order estimates, since we assume a centered dipole 
for the planet's internal field; Nichols et al. [2009] have used auroral observations and a 
more realistic internal field model to show that the northern auroral oval has a radius 
~1.5° smaller than the southern. The model magnetopause shielding field is valid for 
latitudes <50° [Caudal, 1986], and the labeled field lines exceed this limit for radial dis- 
tances <7R S . Increased hot plasma content expands the equatorial distance of the labeled 
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field lines by ~3-4Rs- The ionospheric colatitudes 15—17° are <1° from the observed 
equatorward boundary of the northern auroral oval [Nichols et al, 2009]. Thus inter- 
nal re- configuration of hot plasma, at fixed Rmp, can influence ionosphere-disk magnetic 
mapping and equatorial width of the auroral current layer. A similar effect was noted by 
Grodent et al. [2008] from auroral observations of Jupiter. 

In Fig. 2b, we plot equatorial magnetic and plasma pressure for both models. Field 
strength at p>15Rs for the disturbed disk is more uniform than the quiescent model. 
Total model plasma pressure (disturbed) agrees well, as shown, with that of Sergis et al. 
[2010] (hereafter S10; derived directly from CAPS and MIMI pressures). The radial 
gradient of our profile differs from that of S10, however, which has consequences for 
azimuthal current (see below). Our disturbed disk's maximum plasma j3 of ~6.5 at 
p~14Rs is consistent with SlO's /3 m a:r~3-10 for p~ll-14Rs. 

Fig. 2c compares equatorial volume forces. Again, the smaller gradient in the disturbed 
disk's magnetic pressure for p>15Rs is evident. The corresponding increase in total 
plasma pressure gradient is mainly balanced by increased curvature force of the more 
disk-like field. In the outer model disk, plasma pressure gradient becomes smaller than 
centrifugal force within a few R s of the magnetopause. This result, arising from our im- 
proved pressure representation, differs greatly from AlO's "transition distance" of ~12 R s , 
but agrees with Arridge et al. [2007], who deduced plasma pressure from Cassini magnetic 
data during current sheet encounters. 

Fig. 2d compares the modeled equatorial azimuthal current density, associated with 
pressure gradient and centrifugal force, with SlO's data-derived values. Considering cen- 
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trifugal (inertial) current, the model is in reasonable agreement with the data for p>6Rs. 
The discrepancy is mainly due to the difference between our approach (see A 10) and 
S10 for computing currents. For centrifugal current, we assume AlO's polynomial fit to 
plasma angular velocity co obtained from high-energy particle anisotropies [Kane et al, 
2008], combined with inner magnetospheric values from Wilson et al. [2008]. Combining 
our magnetic field model and profile for ion flux tube content produces the centrifugal 
term. S10 employed CAPS density and velocity moments from Wilson et al. [2008] and 
McAndrews et al. [2009], and their own CAPS ion moment computations, in order to 
derive and fit centrifugal volume force as a function of p. The ratio of our centrifugal 
force to that of S10 ranges from ~2 at ~6Rs to ~0.5 at ~12Rs. 

Considering pressure gradient, we recall (discussion of Fig. 2b) considerable differences 
in plasma pressure gradients between our model and SlO's values. Thus for p>8Rs, our 
disturbed model's total pressure current seems reasonably consistent with SlO's lower 
limits. Current density exceeds inertial at ~8Rs in both model and data. For p<8Rs, 
Sergis et al. [2010] find >25% of total current arises from pressure anisotropy. This is not 
included in our model at present. 

In Fig. 3a,b, we compare model outputs with hot plasma and magnetic pressures for 
CassinVs Rev 40 orbit from the MIMI and Cassini magnetometer (MAG) [Dougherty 
et al., 2004]. The model assumes constant hot pressure Ph along field lines, hence the 
change in local pressure with altitude Z (from the equator) is equivalent to the change 
in the value Pho at the equatorial crossings of the field lines sampled by the spacecraft. 
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Values of Rmp were determined by A10, who modeled individual field components during 
this orbit. 

Hot pressures near the nightside current sheet crossing (CSX1) are well captured by a 
quiescent disk model, while the dayside crossing (CSX2) requires an "intermediate" state. 
Sharp "dropouts" in pressure near days ~67.5 and ~70 (decimal day of year) are also seen 
in the model, and coincide with the spacecraft's exit from and entry into the polar cap, 
characterized by open magnetic flux and relative absence of plasma. This region manifests 
in the model as field lines which do not cross the equator inside the magnetopause, hence 
are assumed to be devoid of plasma. Polar cap boundary crossings are also seen in 
nightside models, near days ~63.9 and ~65.9. For these local times, however, observed 
pressures are small, near noise level, and could obscure any dropout signature. Quasi- 
periodic "pulsations" are present in plasma and magnetic data, which we emphasize with 
vertical lines spaced by the nominal 10.75 hr planetary rotation period. These features are 
due to the quasi-periodic "camshaft" oscillations in the Kronian magnetosphere, whose 
period agrees with that of the Saturn kilometric radiation (e.g. Southwood and Kivelson 
[2007]; Kurth et al. [2007]; Provan et al. [2009]). 

The change in magnetic pressure when crossing the current sheet is greater at CSX1 
than CSX2. CSX1 occurs at p~15Rg in the "middle magnetosphere" where the field lines 
are most disk-like, hence produce a strong contrast in field strength between the inside and 
outside of the sheet (Fig. 2a). CSX2 occurs at p~20Rs towards the "outer region" where 
field lines are more "blunt" (larger radius of curvature), this contrast is thus reduced, and 
the cold plasma sheet is thicker (A 10). Models with more disturbed disk and / or larger 
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Rmp show a change in magnetic pressure comparable to that seen at CSX1. However, the 
observed sheet transit is more rapid than the model analog, perhaps due to plasma sheet 
motion linked to camshaft oscillations. 

Fig. 3c is a CAPS electron spectrogram for the same orbit. The decrease in cold electron 
plasma (counts) with distance from the equator is more pronounced near CSX1 than 
CSX2. This aspect reflects global force balance, which allows the cold population to 
maintain higher pressure for a given Z in regions where the field lines have larger curvature 
radii, and centrifugal confinement is weaker (e.g. A 10, Caudal [1986]). Nightside data are 
better matched by a more quiescent disk and larger Rmp, both possible signatures of a 
thinner nightside plasma sheet, as described by Sergis et al. [2009]. The magnetosheath 
population seen after the outbound magnetopause crossing (MPX, early on day 70) shows 
peak energies similar to those near the exit from the polar cap (via the dayside cusp) at 
day ~67.5. Trajectory information in Fig. 3d includes northern ionospheric colatitude 6>, 
magnetically conjugate to CassinVs location. Regions with ^<13° are shown as a solid 
curve (also shown in Fig. 3c), and represent a model polar cap, located on field lines 
poleward of the auroral oval boundary [B adman et a/., 2006]. 

4. Conclusions 

Saturn's magnetospheric field can be significantly modified by internal changes in hot 
plasma pressure. Magnetic mapping between ionosphere and disk is also dependent on 
internal plasma state (as well as external solar wind pressure, e.g Bunce et al. [2008]). 
This result suggests a corresponding influence on the distribution of currents which flow 
between these regions and their associated auroral emissions. Comparing our models with 
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both field and plasma data allows identification of ring current state and its magnetic 
signatures. 
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Cassini MIMI / CAPS Equatorial Data 




Figure 1. (a), (b) Perpendicular and parallel temperature moments (with respect 
to assumed-southward magnetic field) for protons (H, circles) and water group ions 
(W, crosses), from Wilson et al. [2008] and McAndrews et al. [2009]. Solid curves show 
quartic fits; (c) Hot plasma pressure derived from the MIMI data. Gray squares cover 23 
Cassini orbits, black squares are from Rev 3; (d) Median and quantile pressures computed 
for bins in (e). Solid curves are cubic fits used to calculate PhY a (see text); (e) Radial 
distance bins and number of samples used. 
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Saturn Magnetodisc Models, R MP =25 R g 
Saturn Disc: Quiescent RC Saturn Disc: Disturbed RC 




Figure 2. Model outputs for magnetopause radius Rmp = 25 Rs, and for quiescent 
(left) and disturbed (right) ring current state, (a) Magnetic field lines, partially labeled 
with the colatitude 9i (in degrees) of their northern ionospheric footpoint (2° steps in #j); 
(b) Equatorial pressure. Color code indicates magnetic or plasma pressure (cold, hot and 
total), dashed black curve is from S10; (c) Equatorial volume forces, colored according 
to physical origin. Line style indicates radially outward (solid) or inward (dashed) force; 
(d) Equatorial azimuthal current density, colored according to physical origin, compared 
with the medians (symbols), 10% and 90% quantiles (vertical bars) of the data-derived 
currents of S10 (1 Rs window). 
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Disc Models Compared with Cassini Plasma / Field Data 
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Figure 3. (a) Hot plasma pressure versus time for Cassini 's Rev 40 orbit (symbols are 
medians, vertical bars show 25% and 75% quantiles, for a one- hour window) . Disk model 
pressures are overlaid (see code in panel (b)); (b) Observed and modeled magnetic pressure 
for Rev 40, with MAG data shown as symbols and model curves color- and line-coded 
according to magnetopause radius and ring current activity. Thin vertical dashed lines 
are spaced by the planetary rotation period; (c) CAPS electron spectrogram (one-hour 
average from Anode 5), with current sheet (CSX) and magnetopause (MPX) crossings 
marked. Horizontal line shows the model polar cap interval; (d) Spacecraft trajectory, 
showing cylindrical and spherical radial distance p and r, vertical coordinate Z, local 
time SLT and model colatitude 9i of northern magnetic footpoint. 
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